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ABSTRACT 



We use a purely data-driven rectified factor analysis to identify early-type galaxies with 
recent star formation in Data Release 4 of the Sloan Digital Sky Survey Spectroscopic Cat- 
alogue. We compare the spectra and environment of these galaxies with those of 'normal' 
early-type galaxies, and a sample of independently selected E+A galaxies. We calculate the 
projected local galaxy surface density from the nearest five and ten neighbours (£5 and £10) 
for each galaxy in our sample, and find that the dependence on projected local density, of 
qq ■ the properties of E+A galaxies, is not significantly different from that of early-type galaxies 

C^) with young stellar populations, dropping off rapidly towards denser environments, and flat- 

\Q , tening off at densities 5 0. 1 —0.3 Mpc~ 2 . The dearth of E+A galaxies in dense environments 

■ confirms that E+A galaxies are most likely the products of galaxy-galaxy merging or interac- 

tions, rather than star-forming galaxies whose star formation has been quenched by processes 
unique to dense environments, such as ram-pressure stripping or galaxy harassment. We see 
a tentative peak in the number of E+A galaxies at £10 ~ 0.1 — 0.3 Mpc~ 2 , which may 
represent the local galaxy density at which the rate of galaxy-galaxy merging or interaction 
rate peaks. Analysis of the spectra of our early-type galaxies with young stellar populations 
suggests that they have a stellar component dominated by F stars, ~ 1 —4 Gyr old, together 
with a mature, metal-rich population characteristic of 'typical' early-type galaxies. The young 
stars represent £10% of the stellar mass in these galaxies. This, together with the similarity 
of the environments in which this 'E+F' population and the E+A galaxy sample are found, 
suggests that E+F galaxies used to be E+A galaxies, but have evolved by a further ~ one to 
a few Gyr. Our rectified factor analysis is sensitive enough to identify this hidden population, 
which allows us to study the global and intrinsic properties of early-type galaxies created in 
major mergers or interactions, and compare them with those early-types which have had the 
bulk of their stars in place since a much earlier epoch. 

Key words: Galaxies: evolution; Galaxies: elliptical and lenticular; Galaxies: stellar content; 
Galaxies: starburst; Cosmology: observations. 



1 INTRODUCTION 

There is still much uncertainty about the relative epochs of mass 
assembly and star formation in early-type galaxies (e.g. de Lucia et 
al. 2006; Cimatti, Daddi & Renzini 2006). The most massive galax- 
ies contain the oldest stellar populations, but hierarchical structure 
formation requires that their present-day mass was assembled at a 
recent epoch. On the other hand, galaxy 'downsizing' (e.g. Cowie 
et al. 1996; Treu et al. 2005) means that less massive galaxies have 
more recent star formation. 

One would expect that different formation histories would 
give rise to different observed characteristics in early-type galax- 
ies. Several studies (e.g. Faber et al. 1997; Emsellem et al. 2004; 
Khochfar & Burkert 2005) have in fact found that various sub- 



classes of early-type galaxy exist (e.g. 'core' and 'cuspy' galax- 
ies, fast and slow rotators). We would like to know how these sub- 
classes are related to the various modes of formation of early-type 
galaxies, and to the effects of their immediate environments. 

The stellar populations of early-type galaxies hold the key to 
understanding the difference between early-types created via dif- 
ferent evolutionary routes. If we can identify those galaxies which 
have recently undergone a strong starburst, and we associate that 
recent starburst with a major galaxy-galaxy merger or interaction, 
then we can identify the products of recent merger or interaction, 
map their evolution with environment and redshift, and study their 
intrinsic properties in comparison with more massive 'primordial' 
early-types. 

In this work, we are particularly interested in the relationship 
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between E+A galaxies and 'typical' early-types. E+A galaxies are 
defined as having strong hydrogen Balmer absorption lines, but 
no corresponding [Oil] or Ha emission lines. Their spectra, there- 
fore, resemble the linear superposition of a characteristic early-type 
galaxy spectrum and that of an A star. This suggests that these 
galaxies have undergone a burst of star formation in the last Gyr 
or so, but that there is no on-going star formation. It is hypoth- 
esised that these galaxies will eventually settle down as 'normal' 
passively-evolving early-types. 

E+A galaxies may therefore represent an intermediate stage 
between gas-rich, disc-dominated, rotationally-supported, star- 
forming systems on the one hand, and bulge-dominated, pressure- 
supported, passive systems on the other. Hence, understanding the 
physical processes by which the star formation is triggered and sub- 
sequently truncated, is a vital stage in understanding the evolution 
of early-type galaxies. 

There are three possible scenarios which could explain the 
existence of E+A galaxies. The first is that they are in fact dust- 
obscured star-forming galaxies. This can be verified by observing 
at radio wavelengths, where star formation can be detected, but 
dust obscuration is not significant. Miller & Owen (2001) and Goto 
(2004) observed between them a total of 51 E+A galaxies. Moder- 
ate star formation was detected in only two of these galaxies, so it 
is unlikely that the majority of E+As are dusty starbursts. 

More plausibly, then, E+A galaxies may be the product of 
disk-disk merging or interaction, at a stage when the starburst 
has burnt itself out, or they may be galaxies whose star forma- 
tion has been initiated and then quenched by, for example, ram- 
pressure stripping or tidal stripping. The environment of these 
galaxies therefore holds the key to understanding them. Although, 
at intermediate redshift (0.33 < z < 0.83), Tran et al. (2003) find 
that E+As make up a non-negligible fraction (~7-13%) of cluster 
members, suggesting a cluster-related mechanism for the truncation 
of star formation, such as ram-pressure or tidal stripping, at lower 
redshift (z ~ 0.1), E+As appear to occur in lower-density environ- 
ments, indicative of merging or interaction processes (Zabludoff 
et al. 1996; Blake et al. 2004). In addition, both Zabludoff et al. 
(1996) and Blake et al. (2004) found a significant number of their 
E+A galaxies displayed tidal features and disturbed morphologies, 
typical of merger remnants. More recently, Goto (2005) found that 
E+As selected from the Sloan Digital Sky Survey (SDSS) Data 
Release 2 have an excess of local galaxy density only at scales of 
< 100 kpc, and not at the scale of either clusters or larger-scale 
structures. This is also highly suggestive of a merger or interaction 
origin for E+As. 

Conventional methods of identifying E+A galaxies use mea- 
surements of the equivalent width of hydrogen Balmer absorption 
lines (typically Ha and US) and [Oil] emission. However, this can 
suffer from dust obscuration of the emission lines, and emission- 
filling of H Balmer lines, and it requires a high signal-to-noise (> 
10 per pixel). In addition, there has been some debate over the cri- 
teria for selecting E+A galaxies, with some authors choosing Ha, 
H<5 and [Oil], with differing line strength criteria (e.g. Goto 2004, 
2005), and some including H7 as well as Ha, US and [Oil] (e.g. 
Zabludoff et al. 1996; Blake et al. 2004). 

Here, we use a purely data-driven rectified factor analysis data 
model, developed in Nolan et al. (2006), on the spectra of early- 
type galaxies in the SDSS Data Release 4 (DR4). This was specif- 
ically designed to exploit large data sets in a model-independent 
way. It uses all the data present in the spectra, i.e. the spectral fea- 
tures and the overall shape of the continuum, to quantify the rel- 
ative strength of recent star formation in early-type galaxies ex- 



tracted from the database, in a rapid, objective and robust manner. 
Our method allows us to reconstruct even partial spectra, without 
the need to measure individual lines accurately. The advent of large 
surveys, such as the Sloan Digital Sky Survey allows us to study 
large samples of these rare objects (< 1 % of the overall zero- 
redshift galaxy population, Zabludoff et al., 1996), in a way that 
has not been previously possible. 

It has recently been shown that E+A galaxies occur more fre- 
quently at high redshift (1 < z < 2) compared with the local uni- 
verse, at least amongst the massive red galaxy population (Doherty 
et al. 2005, Le Borgne et al. 2006). They may, therefore, represent 
a critical phase in the evolution of early-type galaxies. However, it 
is difficult to constrain the properties of these objects at high red- 
shift, due to the difficulty of measuring absorption lines accurately 
at that redshift, and in obtaining statistically significant spectro- 
scopic samples. Studying large samples at low redshift is therefore 
vital for the understanding of these objects. 

In this work, we use our data modelling to identify early-type 
galaxies (ETGs) which have a young stellar population, suggesting 
that they have undergone a recent burst of star formation. We com- 
pare these galaxies with E+A galaxies, selected in a completely 
independent manner, and investigate the environment of our galaxy 
samples, in order to understand the physical mechanisms leading 
to the cessation of star formation. 

In Section [2] we define the selection of the early-type sam- 
ple, and the E+A subsample for comparison. We discuss the factor 
analysis of the early-type spectra that is used to reveal young stel- 
lar populations in Section[3] In Section[4] we compare our results 
from the analysis of early-type galaxies with the sample of E+A 
galaxies. In Section [5] we investigate the relationship between the 
galaxies in our various samples and their local environment. Our 
conclusions are discussed in Section[6] 

We use H = 70 kms^Mpc" 1 , ft m =0.3 and f2 A =0.7 
throughout. 



2 THE SAMPLE 

Our sample of early-type galaxy spectra was selected from the 
SDSS DR4 spectroscopic catalogue. We follow the criteria of 
Bemardi et al. (2003) to select early-type galaxies: 

• the likelihood that a de Vaucouleurs profile fits the radial sur- 
face brightness distribution is at least 1 .03 times the likelihood that 
the radial profile is exponential; 

• the concentration index rgo /rso > 2.5 in the i-band; 

• the r-band Petrosian magnitude < 17.5 

• the photometric signal-to-noise > 4 

• the galaxy is not blended, de-blended or a child galaxy 

• a surface brightness profile has been fitted, and it is not satu- 
rated 

• redshift, z < 0.3 and the confidence in the redshift, zConj> 
0.95. 

We also select E+A galaxies from the SDSS DR4. Here, we 
use the criteria of Goto (2005), so that we can compare the results 
of our factor analysis of the early-type galaxy sample with their 
sample of E+As. The criteria are: 

• Ha equivalent width (EW) < 3.0 A 

• US EW < - 5.0 A 

• [OII]EW<2.5A 

• S/N > 10 
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• z > 0.01, 

where emission lines have positive values. The sample is not re- 
stricted by any morphological criteria. 

This gives us a sample of over 26,000 early-type galaxy 
spectra, and 765 Goto E+As, which we de-redshift and re-bin at 
the SDSS full-width half-maximum (FWHM) spectral resolution, 
with the observational errors supplied re-binned correspondingly 
in quadrature. 

In this work, we consider only those galaxies in our samples 
which lie in the range 0.06 < z < 0. 14. Around 50% of the galaxies 
in our ETG sample lie within this range, and restricting ourselves 
to a limited redshift range minimises any effects arising from the 
fixed size of the SDSS fibre diameters. At z = 0. 1, the central ~ 
5.5 kpc are covered, which represents 20—40% of the total stel- 
lar light, which should be a good representation of the bulk stellar 
light. The emission lines used to select E+A galaxy spectra are 
measured from the same spectra that we decompose using our data 
model, so we are of course comparing results extracted from ex- 
actly the same region of each galaxy. 

The redshift cut means that, in this work, we do not explore 
the evolution of ETGs with young stellar populations with redshift. 
However, we also need not worry about fc-corrections to the magni- 
tudes. The redshift range 0.06 < z < 0. 14 is less than 1 Gyr of stel- 
lar evolution; fc-corrections would probably introduce uncertainties 
at least on this scale. Hence, we use a magnitude limit of r < 17.5 
for all our samples, where r is the Petrosian r-band magnitude. The 
SDSS DR4 spectroscopic catalogue is complete to this limit. 



3 MODELLING THE SPECTRA 

In order to analyse the spectra of our early-type galaxy sample, 
we model the data using the rectified factor analysis data model 
that we developed in Nolan et al. (2006). This model-independent 
technique is able to recover physically meaningful components of 
the observed spectra in a timely manner. Our approach and the un- 
derlying rationale for choosing the rectified factor model for mod- 
elling the data are described in detail in Nolan et al. (2006, §4.1), 
along with the details of an implementation employing variational 
Bayesian techniques (Nolan et al. 2006, §3.1). This modelling is 
an unsupervised technique, i.e. a data-driven exploratory analysis 
tool, designed to help us discover the underlying structure of the 
data and to derive features for either interpretation or subsequent 
quantitive analyses. Unsupervised techniques are machine learning 
approaches, which do not require a-priori specified targets for the 
training examples but instead, the data density is modelled. This 
model also includes latent variables that are then used to infer pat- 
terns in the data that are hoped to reveal useful relationships. 

The model is first trained on half of the SDSS early-type 
galaxy spectra (~ 13000 galaxies), and the recovered components 
are then used to reconstruct the remaining ETG spectra, and the 
spectra in the Goto E+A sample (GE+A). Following Nolan et al. 
(2006, see §4 and Fig. 4), we chose to have two spectral compo- 
nents, which was shown to be sufficent to represent the bulk of the 
stars in early- type galaxies. Fig. Q] shows the two component spec- 
tra recovered from the factor analysis, together with a 10 Gyr, 2.5 
Zq single stellar population model (Bruzual & Chariot 2003) and a 
typical super-solar F star spectrum (Santos et al. 1995) for compari- 
son. Features typical of a mature, metal-rich stellar population, and 
a young stellar population are marked on the plot. It can be seen 
that the first component represents a contribution dominated by 
younger stars (~l^t Gyr for a typical F star at the main sequence 
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Figure 1. From top to bottom: a 10 Gyr, 2.5 Zq single stellar population 
model (Bruzual & Chariot 2003, black); component 2 from our analysis 
(red); a super-solar metallicity F star spectrum (Santos et al. 1995 green); 
component 1 from our analysis (blue). The dotted lines mark some of the 
absorption features in the spectra which are typically strong in young stel- 
lar populations, and the dashed lines mark some of the absorption features 
which are typically strong in old, metal-rich stellar populations. From left 
to right, the absorption line species are: He (3970 A), H<5 (4102A), H7 
(4340A), Bf3 (4861 A), Mgb (5175A), NaD (5893A), Ho, TiO (7126A). It 
can clearly be seen that the first component (blue) represents a younger stel- 
lar contribution (green), whereas the second component (red) represents a 
mature stellar population (black). 

turn-off point). It has strong hydrogen B aimer absorption lines and 
a corresponding lack of evolved metal absorption lines. The sec- 
ond component has the characteristic shape and spectral features of 
a mature, metal-rich stellar population, such as those that are ex- 
pected to dominate in early-type galaxies. This component has a 
well-developed 4000A break, typical of such populations. Clearly, 
the two components recovered represent young and old stellar con- 
tributions to the spectra of ETGs, although of course they cannot be 
exactly interpreted as single stellar populations. Although the sec- 
ond component is well-fit by the 10 Gyr stellar population, the first 
component is slightly redder than that of the F star. This probably 
reflects both a more complex star formation history than a simple 
two-component model, and that the 'real' spectrum is that of the in- 
tegrated young stellar population, rather than that of a single stellar 
type. Ferreras et al. (2006), using a principal components analysis 
on the 3500— 7500A spectra of 30 ellipticals in groups and the field, 
recover similar representations of old and young stellar populations 
in their first two components. Their results are also consistent with 
the representation of the bulk of stars in early-type galaxies by two 
components. 



4 COMPARISON OF THE EARLY-TYPE GALAXIES 
AND THE E+A SAMPLE 

Within our chosen redshift range, we are left with 13,111 ETGs. 
Keeping only those galaxies in the GE+A sample which have Ha, 
H/3 and [Oil] securely determined leaves us with 196 GE+As. 
Within the ETG sample, we define a sub-sample of ETGs with 
young stellar populations (YSPETs). These are galaxies for which 
al, the contribution (or weight) of the 'young' data model compo- 
nent is greater than a2, the weight of the 'old' data model compo- 
nent. By comparison with stellar population models, we estimate 
al > a2 is equivalent to an approximate lower limit to the young 
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stellar population of £10% by stellar mass. YSPETs account for 
17% of the ETGS (2245 galaxies). 

In Fig. |2] we plot the mean spectra for each of these three 
catagories, with the vertical normalisation adjusted for ease of com- 
parison. The mean GE+A and YSPET spectra are very similar 
longwards of H<5, and both have stronger hydrogen Balmer absorp- 
tion lines and a steeper 5000— 8000A slope than the ETG spectrum. 
This is characteristic of the presence of contributions to the flux 
from young stars. The GE+A mean spectrum shows contributions 
from the characteristic features of even younger stars than does the 
YSPET mean spectrum. It has stronger hydrogen Balmer absorp- 
tion lines, enhanced flux just longwards of 4000A, and the flux falls 
off more steeply towards the blue end of the spectrum shortwards 
of 4000 A (see Fig. [U. 

In Fig. [3] we show the likely stellar composition of the GE+A 
and YSPET galaxies. The E+A galaxies, as one would expect, re- 
semble the linear combination of a mature, metal-rich stellar popu- 
lation (the 'E' component) together with a population dominated by 
young (~ 0.7 Gyr) A stars. The mean YSPET spectrum is similar, 
but in this case, the younger stellar component is better-represented 
by F stars, which dominate the flux of a young stellar population at 
a later epoch than A stars (~ 1— 4 Gyr). Hence, our factor anal- 
ysis has allowed us to define a sample of early-type galaxies with 
young-intermediate aged stars, which, as F stars dominate the in- 
tegrated spectrum of a stellar population for longer than the short- 
lived A stars, presents us with a larger sample of galaxies which 
have undergone recent star-forming activity than E+A samples are 
able to. 

Furthermore, in Fig. [4] we plot the number densities of the 
weights of the data model components for our three classes of spec- 
tra. We have further divided the GE+A galaxies into those that 
would be classified in our scheme as early-types, and those that 
would not. Only ~ 17 % of the GE+As (34) are typical early-type 
galaxies, of which 2/3 are in the al > a2 region. The distribution 
of the GE+As is less peaked in either al or a2 than the YSPETs, 
which broadly follow the same distribution as the ETGs. 

That most of the GE+A galaxies are not ETGs is consistent 
with the idea that they are the result of galaxy-galaxy mergers or 
interactions, and are at an earlier epoch in their post-merger evolu- 
tion than the YSPETs. Their young stellar component is consistent 
with being less than 1 Gyr old. This is less than the typical time 
for the colours and surface brightness distributions of merger rem- 
nants to relax into those which are typical of ETGs, and which we 
use to identify ETGs in the SDSS catalogue. Of course, it should 
also be noted that if GE+As are quenched spirals, whose star for- 
mation has been recently halted by ram-pressure or tidal stripping, 
one would also expect them not be classified as ETGs under our 
scheme. 



5 THE ENVIRONMENTAL DEPENDENCE 

In order to investigate the mechanisms via which E+As are formed, 
we have quantified the environment of our sample galaxies. First, 
we find CT10 and 0-5, the projected distance between the galaxy and 
its 10th and 5th nearest neighbours in the SDSS spectroscopic cat- 
alogue. We use a Petrosian magnitude limit, r < 17.5, and do not 
fc-correct our magnitudes, due to the limited redshift range of our 
sample (see S|2}. The SDSS DR4 spectroscopic survey is ~ 92 % 
complete (Ani Thakar, private communication), so we are able to 
take advantage of the spectroscopic redshifts in determining envi- 
ronment. 
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Figure 2. Mean spectra for the whole ETG sample (black, bottom), the 
YSPET sub-sample (red, middle)), and the GE+A sample (blue, top). The 

spectra are normalised to unity at 5500 A and the vertical normalisation 
adjusted for comparison. Longwards of H<5, the E+A and the YSPET pop- 
ulations look very similar. The mean ETG spectrum is flatter, and does not 
have such pronounced hydrogen Balmer absorption lines. The dotted and 
dashed lines are the same as for Fig.fT] 
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Figure 3. LHS, from top to bottom: A 10 Gyr, 2.5 Zq stellar population 
model (Bruzual & Chariot 2003, red:); An "A star" spectrum (Kurucz 1993, 
blue); an E+A spectrum, which is a linear combination of the first two 
spectra (green); and the mean GE+A spectrum (black). RHS: as for the 
LHS, but here we use a super-solar metallicity "F star" (Santos et al. 1995, 
blue) instead of the A star, and the lowest (black) spectrum is the mean 
spectrum of the YSPET population. The mean GE+A spectrum, as one 
would expect, is very similar to our constructed E+A spectrum. The mean 
YSPET spectrum, however, looks more like a 'E+F' spectrum. The dotted 
and dashed lines are as for Fig. □ 

We cut in peculiar velocity at ± 1000 kms~ from each 
galaxy, to attempt to avoid background or foreground objects, but 
include cluster members if the target galaxy belongs to a cluster. 
We have not attempted here to determine cluster membership, and 
a more sophisticated analysis is reserved for a later paper. Here, we 
are only interested in the projected local galaxy density. 

The projected distances are then transformed into an estimate 
of the local projected surface density, £10 and E 5 . We have made 
a cut in a n to allow for the effects of the survey boundary, which 
could lead to an underestimate of local density for those galaxies in 
low-density regions. Following Clemens et al. (2006) and Miller et 
al. (2003) we cut at a n (max) = 7.2 Mpc. 

In Fig. [5] the normalised number distributions of local galaxy 
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Figure 4. Distribution of the weights al (the weight of the 'young' component from our data modelling) and a2 (the weight of the 'old' component from our 
data modelling), for: all the early-type galaxies in the sample (black); Goto E+As (blue, cross-hatched ); YSPETs (red, hatched). In the left plot, the GE+A 
sample is further divided into those that would be classified in our scheme as early-types (green) and those that would not (yellow). 



surface density, Eio and Eio, are shown for the YSPET and the 
GE+A galaxy populations. The similarities between the samples 
can be clearly seen. The distribution of the GE+As follows that 
of the YSPETs in both cases. The fraction of GE+As perhaps 
peaks at around Eio ~ 0.1 — 0.3 Mpc" 2 , and both populations fall 
off strongly towards higher density environments. The fraction of 
YSPETs flattens off after this peak in Eio, and in E 5 , the YSPET 
distribution is consistent with a flattening off at surface densities 
less than 0.2-0.3. 

Performing a Kolmogorov-Smirnov test on the GE+A and 
YSPET populations returns the probabilities, (p — 1) = 4 x 10 -4 
(E5) and 0.3 (Eio) that the two samples are taken from different 
distributions. Therefore, in both cases, despite our completely in- 
dependent selection procedures, the distributions of the two popu- 
lations are not significantly different, consistent with the hypothesis 
that our YSPETs are in fact evolved E+A galaxies ('E+F' galax- 
ies). The decreasing fraction of E+As and YSPETs towards in- 
creasing densities strongly argues for the hypothesis that the ma- 
jority of GE+As or YSPETs are not the product of cluster-related 
processes (e.g. ram-pressure stripping). This confirms the result 
of Goto (2005) that E+A galaxies are likely to be the product of 
galaxy-galaxy mergers and interactions in a group environment. 

Studies of star-forming galaxies of all morphologies in the 
SDSS and the Two Degree Field Galaxy Redshift Survey (2dF- 
GRS) (e.g. Lewis et al. 2002; Gomez et al. 2003; Balogh et al. 
2004; Poggianti et al. 2006) find that the star formation rate falls 
off steeply with increasing density, but converges to the field galaxy 
distribution at a projected density of ~ 1 Mpc" 2 . In addition, 
Balogh et al. (2004) find that the fraction of star-forming galaxies 
of all morphologies falls steeply with respect to E5 ~0.4 Mpc" 2 , 
flattening off at densities less than this, in good agreement with our 
results (see their Figs 5 and 6). 

The peak in fraction at Eio ~ 0.1 — 0.3 Mpc" 2 may indicate 
the optimum local galaxy surface density for galaxy-galaxy merg- 
ers and interactions, representing a trade-off between sufficient 
density for frequent galaxy-galaxy encounters, and slow enough 
relative velocities that those encounters lead to an interaction. In- 
vestigation of the results of N-body simulations would be able con- 
firm or reject this. However, there is support for this argument in the 
recent results of Ferreras et al. (2006). They find a higher mass frac- 
tion of younger stars in some group galaxies than in the field galax- 
ies, implying a more complex star formation history for galaxies in 
groups than in the field, in qualitative agreement with our tentative 
peak of recent starburst activity. 

Further evidence for a peak in star-formation with respect to 
local galaxy density is found by Porter and Raychaudhury (2006, 



in preparation). Along filaments joining rich clusters, they find a 
peak in star formation at 3—4 Mpc from cluster centres. However, 
a careful assessment of group, cluster and filament membership is 
necessary, as the precise location of the peak depends upon whether 
a galaxy is a member of any these environmental classes. 

Many authors have found that the star-formation rate begins 
to fall at densities Zl Mpc -2 (e.g. Lewis et al. 2002; Gomez et al. 
2004; Balogh et al. 2004; Poggianti et al. 2006; Porter and Ray- 
chaudhury 2006, in preparation). However we find, as do Balogh et 
al. (2004), that the fraction of star-forming galaxies starts to fall off 
at lower densities (0.1—0.4 Mpc" 2 ), with the caveat that peaks and 
breaks with respect to density depend on a measurement of envi- 
ronment more complex than simple projected local galaxy density, 
as mentioned above. This is suggestive that actively star-forming 
disk galaxies have their (moderate compared with disc-disc merg- 
ing starbursts) star-formation quenched in group environments, de- 
creasing the active-to-passive galaxy ratio. At the same time, more 
disc-disc merging or interaction occurs in the group environment, 
with associated massive star formation, keeping the star-formation 
rate high at higher densities, whilst contributing to the fall in the 
number of disc galaxies. Hence, both morphological and spectro- 
scopic data are required in order to study these different processes 
effectively. 



6 CONCLUSIONS 

We have successfully used our data-driven rectified factor analy- 
sis to identify early-type galaxies with young stellar populations, 
without having to apply the far more computationally expensive 
method of fitting detailed stellar population synthesis models to 
a sample of > 13,000 ETGs in the redshift range 0.06 < z < 
0.14. We have used the sample generated as a result to map the 
relationship between early-type galaxies, E+A galaxies and early- 
type galaxies with young stellar populations. Our analysis recovers 
two components from the set of early-type galaxy spectra, one of 
which represents the contribution to the flux from mature, metal- 
rich stars, and the other which represents the contribution from F 
stars, ~1^4 Gyr old. Our method allows us to reconstruct spectra 
even when the observed spectra are incomplete. We do not need 
high signal-to-noise measurements of individual lines. We identify 
those early-types which have undergone recent (~ 1-4 Gyr) star for- 
mation, with the young stars accounting for at least 10% of their 
stellar mass (YSPETs). At redshifts 0.06 < z < 0.14, ~17% of 
ETGs are classified as YSPETs. We propose that these YSPETs 
are 'E+F' galaxies, whose spectra are a linear combination of a 
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Figure 5. Normalised distributions of local surface density: early-types with a2 > al (red, solid) and Goto's E+As (blue, cross-hatched) and £5 (top), £10 
(bottom). The populations are normalised so that the area under each line represents 100% of the galaxies in that sample. 



typical early-type galaxy spectrum and that of an F star, analogous 
to E+A spectra. 

We also select a sample of E+A galaxies in a completely in- 
dependent way, without using the morphological criteria we used 
to identify early-type galaxies, but instead relying on the emis- 
sion and absorption-line criteria of Goto (2005). The characteristic 
'E+A' spectra of the Goto-selected sample, and the fact that many 
of these E+As do not fulfil the colour and concentration criteria 
for typical early-type galaxies, suggests that these galaxies are in 
an early (< 1 Gyr) post galaxy-galaxy merger or interaction phase. 
Goto (2005) notes that most of the E+As in his sample have con- 
centrations as high as typical ellipticals {Cm < 0.4), so it is likely 
that many of these galaxies are not included in our early-type sam- 
ple because they are too blue, indicating that massive (A star and 
earlier) stars are still burning. If E+As are recently quenched star- 
forming galaxies, we would expect them to be located preferen- 
tially in dense cluster regions, whereas we see the majority resid- 
ing in field or group environments, and the number falls off towards 
denser environments, in agreement with Goto's result. Their high 
concentrations and preferentially low-density environments argue 
strongly that these are the products of galaxy-galaxy mergers or 
interactions, which will eventually settle down to become typical 
early-type galaxies. 

We note that there is in fact a slight peak in the local galaxy 
surface density distribution of E+A galaxies at 0.1-0.3 Mpc -2 . 
This may represent the optimum local galaxy surface density for 
galaxy-galaxy mergers and interactions. Realistic cosmological 
simulations should be able to confirm or reject this hypothesis. 

Given that the distribution of local galaxy surface density 
of the GE+A galaxies and that of the independently determined 
YSPET populations are not significantly different, one might ex- 
pect a similar history for both sets of galaxies. We hypothesise 
that the YSPETs, with typical ETG colours and concentrations, 
and 'E+F' spectra, represent an intermediate evolutionary stage 
between E+As and early-types. As this is a longer-lived stage than 
E+A galaxies (F stars persist on the main sequence for longer than 
A stars, and dominate the integrated flux at ages S5 Gyr), we are 
able to identify a larger sample of these galaxies than is possible 
with the very rare E+As. We therefore have identified a powerful 



tool for studying in detail, both globally and individually, the evo- 
lution of early-type galaxies. 

Our future work will compare the dependence on local density 
of E+A and E+F galaxies that we have recovered here, with those 
recovered from simulations. We intend to perform the same analy- 
sis on 2dFGRS galaxies and a sample of high-redshift galaxies for 
comparison. We also intend to investigate the relationship between 
E+A and E+F galaxies by studying their kinematics in detail using 
integral field spectroscopy. 
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